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DIFFUSION  PROBLEMS  IN  BONDED  NONHOMOGENEOUS 
MATERIALS  WITH  AN  INTERFACE  CUT 

by 

Fazil  Erdogan  and  Murat  Ozturk 
Lehigh  University,  Bethlehem,  PA  18015 
ABSTRACT 

In  this  paper  the  mixed  boundary  value  problem  for  a  nonhomogeneous  medium 
bonded  to  a  rigid  subspace  is  considered.  The  main  objective  is  to  investigate  the 
techniques  that  would  lead  to  analytically  tractable  solutions  and  to  provide  examples 
comparing  the  results  of  various  kinds  of  material  nonhomogeneities.  The  problem 
studied  is  a  two  dimensional  diffusion  problem  in  which  the  interface  contains  a  plane 
crack.  An  elastic  medium  under  antiplane  shear  loading  is  used  to  formulate  the 
problem.  However,  the  results  may  be  interpreted  in  terms  of  any  number  of  steady- 
state  diffusion  phenomena.  The  method  used  is  essentially  an  inverse  method  in  the 
sense  that  it  provides  the  material  constitutive  behavior  for  which  the  mixed  boundary 
value  problem  can  be  solved  rather  than  solving  the  problem  for  a  given  material.  Two 
different  methods  are  described  and  some  numerical  examples  are  given. 

1.  Iatnoductioa 

Increasing  concerns  in  recent  years  with  mechanical  failure  initiating  at  the 
interfacial  regions  in  many  technologically  important  multiphase  materials  require  a 
better  understanding  of  the  interaction  between  flaws  that  may  exist  in  these  regions 
and  the  applied  loads  and  other  environmental  factors.  The  conventional  approach  to 
studying  the  thermomechanics  of  such  materials  is  based  on  the  assumption  that  the 
composite  medium  is  piecewise  homogeneous  and  the  flaws  may  be  represented  by  plane 
cuts  or  cracks.  On  the  other  hand  in  most  bonded  materials  the  interfacial  region 
appears  to  have  a  structure  which  is  generally  different  than  that  of  the  adjacent 
materials.  In  many  cases,  such  as  in  plasma  spray  coating,  sputtering,  ion  plating  and 
in  some  diffusion  bonded  materials,  the  thermomechanical  properties  of  the  region  are 
"graded"  in  the  sense  that  the  interfacial  region  is  a  nonhomogeneous  continuum  of 
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finite  thickness  with  very  steep  property  gradients  [lj.  [2j.  Among  other  applications  of 
bonded  nonhomogeneous  materials  one  may  mention  certain  geophysical  materials  with 
naturally  graded  compositions  such  as  shale/sandstone  and  materials  which  have  highly 
temperature  dependent  properties  and  are  under  steep  temperature  gradients.  However, 
perhaps  the  most  important  reason  for  studying  the  diffusion  and  fracture  problems  in 
bonded  nonhomogeneous  materials  is  the  technological  potential  of  the  so-caUed 
functionally  gradient  materials  (FGM).  These  are  the  multiphase  materials,  interlayers 
and  coatings  synthesized  in  such  a  way  that  the  volume  fractions  of  the  constituents  are 
varied  continuously  in  thickness  direction  to  give  a  predetermined  composition  profile 
[3]  -  [9],  The  material  thus  obtained  is  known  to  have  some  wide  ranging  highly 
desirable  properties.  From  the  viewpoint  of  failure  mechanics  some  of  these  properties 
are  reduced  thermal  stresses,  residual  stresses  and  stress  concentration  factors  [3], 
improved  bonding  strength  [G] .  and  improved  toughness  and  corrosion  and  fatigue  crack 
growth  resistance.  Most  of  the  current  research  on  FGMs  appears  to  be  in  the  area  of 
ceramic  coatings  motivated  by  a  variety  of  thermal  shielding  problems  and  in  growing 
superconducting  or  diamond  films  over  homogeneous  substrates  [5]. 

The  general  problem  of  interest  here  is  the  heat  diffusion  and  the  fracture 
mechanics  problems  in  a  medium  that  consists  ot  a  nonhomogeneous  layer  bonded  to  a 
homogeneous  substrate  and  contains  an  interface  crack.  The  primary  objective  is  to 
investigate  the  techniques  which  would  lead  to  analytically  tractable  solutions  of  the 
related  mixed  boundary  value  problems  and  to  provide  some  examples  comparing  the 
results  of  various  kinds  of  material  nonhomogeneities.  To  do  this  we  consider  a  problem 
which  is  relatively  simple  to  work  with  and  yet  has  all  the  analytical  features  that  need 
to  be  investigated.  The  problem  studied  here  is  a  two-dimensional  diffusion  problem  for 
a  nonhomogeneous  medium  bonded  to  a  rigid  half  space  having  a  crack  along  the 
interface.  Even  though  the  terminology  used  will  be  that  of  an  elasticity  problem  for  a 
medium  under  antiplane  shear  loading,  the  results  may  be  interpreted  in  terms  of  any 
number  of  diffusion  phenomena. 

2.  Formulation  of  the  Problem  -  A  Direst  Method 

Consider  the  antiplane  shear  problem  for  a  semi  infinite  nonhomogeneous  elastic 
medium  shown  in  Fig.  1.  It  is  assumed  that  the  medium  is  bonded  to  a  rigid  half  space 
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along  the  y=0  plane,  the  interface  contains  a  crack  on  y=0,  -a<x<a,  and  the  medium  is 
subjected  to  arbitrary  antiplane  shear  loading  which  may  be  mechanical  or  thermal  in 
origin.  It  is  also  assumed  that  the  problem  is  solved  under  actual  loading  conditions  in 
the  absence  of  a  crack  and  through  a  proper  superposition  is  reduced  to  a  local 
perturbation  problem  in  which  the  crack  surface  tractions  are  the  only  non-zero  external 
loads.  Furthermore,  if  we  assume  that  the  shear  modulus  y  is  a  function  of  y  only,  the 
non-vanishing  stress  components  and  the  equilibrium  equation  may  be  expressed  as 


(j  —  u  Sm  u  — 
ry  ;  dx  ’  y:  f  dy  ’ 

r-2  On  pm  „ 

w'TyTy  =  0’ 


(2.1) 

(2.2) 


where  w  is  the  z  component  of  the  displacement  vector.  In  the  crack  problem  under 
consideration  (2.2)  will  have  to  be  solved  under  the  following  mixed  boundary 
conditions: 


ayz  (x,  0)  =  <70  (x)  ,  -  a<  x<  a 
w  (x,  0)  =  0  ,  a  <  jxj  <  oo  , 


cr0  (x)  being  the  known  crack  surface  traction. 
Let  the  solution  of  (2.2)  be  given  by 

u  (x.y)  -  tz  f  F{y,a)  e~'QX  da. 

From  (2.2)  and  (2.4)  it  follows  that 

+M  4^-q2F=0 

dy2  dy 

where 


p(y)  = 


A  y) 
h  (y)' 


If  we  replace  the  unknown  function  F  by  H  so  that 


(2.3) 


(2.4) 


(2.5) 


(2.6) 
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F\y,a)  =  H(y,a)  exp 
equation  (2.5)  becomes 


JJP  U)  dt\  =  H{y,a)  (/i  (y))  2 


(2.7) 


(2.8) 


We  will  now  look  for  a  particular  class  of  functions  p(y)  for  which  (2.8)  has  an 
analytical  solution.  The  simplest  of  such  class  of  functions  are  obtained  by  assuming 
that 


1  j*g  ,  1  -2 

2  dy  '  4  P 


c0  . 


(2.9) 


where  c0  is  a  constant.  Three  classes  of  functions  satisfying  (2.9)  may  thus  be  obtained 
as  follows: 


(a)  c0-  3 2  : 

p{y)  =  =  23  .  p{y)  =  pQcx p  (  =  2 3y)  .  (2.10) 

p{y)  =  23  coth  (3y  +  6)  .  p  {y)  =  p0  sinh2(3y  +  6)  ,  (2.11) 

p(y)  -  23  tanh(3y  4-  <5 )  ,  p(y)  =  pQ  cosh2!Jy  +  6)  ;  (2.12) 

(b)  ca  =  -  32  : 

ply)  =  -23  tan(3y+S )  ,  p  (y)  =  CQSi^y  +  ^  >  (2-13) 

p(y)  =  23  cot(3y  +  8)  ,  p(y)  =  //0  sm2(,dt/-b  6)  ;  (2-14) 

(c)  cQ  =  0  ; 

p(  y)  =  o  -  y(y)  =  Mo  >  (2-15) 

p(y)  =  Jylg  '  MP)  =  ho  4  »2  •  (2-16) 
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where  3,  6,  and  //0  are  arbitrary  constants  (/t0  >  0). 

From  (2.S)  and  (2.9)  it  is  seen  that 

d2  H  \2  ii  __  n  \2 _ 2  ,  „  \  i  2  1/2 

^^2  —  A  rl  —  (J,  A  —  q  +  c0  ,  A  =  (q  t  cQ)  ' 

Formally,  from  (2.4),  (2.11)  and  (2  .17)  we  now  obtain 

,  1  /  a  ( y)  r  x 

w{  £,y)  =  2^V"7Z^“ )  J  "  '4(a)  ex p(  -  Xy  -  iax )  da  , 


(2.17) 


(2.13) 


where  the  constant  yQ  is  introduced  for  dimensional  considerations  and  the  function  .4  is 
unknown  to  be  determined  from  (2.3).  The  physics  of  the  problem  requires  that  w  -*  0 
for  y  -*  oc.  Thus,  from  (2.17)  and  (2.13)  it  may  be  observed  that  A  defined  by 
A  =  (a2  +  c0)1/,:2  must  be  positive  for  all  values  of  a  and,  therefore,  the  negative  values  of 
c0  or  shear  moduli  given  by  (2.13)  and  (2.14)  would  not  be  acceptable  for  the  solution 
expressed  by  (2. IS).  From  (2.1)  and  (2.13)  the  stress  components  are  found  to  be 


y.  ( y.  ( y  A  1  / '  f  x 

°zz  =  V^ThT J  j  -  x  iaA[u)  cxp(  ~  Xy~  lQX)  da  > 


*  =  ~  pt  (jT(y))  J  -  3o  (A/1  +  1  M')  eXp(  ~Xy~  *0fI)  da- 


(2.19) 


(2.20) 


By  substituting  from  (2. IS)  and  (2.20)  into  (2.3)  one  would  obtain  a  pair  of  dual 
integral  equations  to  determine  .4(a).  The  problem  may  also  be  reduced  to  an  integral 
equation  by  defining 


q~x  =  g{x) 


(2.21) 


From  (2.3b),  (2. IS)  and  (2.21)  it  may  thus  be  seen  that 


ia.4  (a)  =  —  f  g(t )  e,c 

J  -  a 


(2.22) 


Substituting  now  from  (2.20)  and  (2.22)  into  (2.3a)  we  obtain 


lim  J 
y~*  +  02 


^  J  <?(t)dt  J  K(y.a)eiai't  ~  da  =  (7yz(x,0)  =  a0(x),  —a<x<a 


(2.23) 


where 


K(y,a)  = 


My)  +  Ay)  A  _A„ 

- - e  j 

ia 


(2-24) 


To  investigate  and  to  separate  a  possible  singular  part  of  the  kernel  in  (2.23),  the 
asymptotic  behavior  of  the  inner  integral  for  j  a  j  -+  oo  must  be  examined.  From  (2.17) 
it  is  seen  that  A  -+  |  a  |  for  |a|  -*  ec.  Thus,  for  ja|  -»  oo  from  (2.24)  we  obtain 


i*(y)  la 

i  * 


e 


(2.25) 


By  adding  and  subtracting  Kx  in  (2.23)  and  by  observing  that 


00  1  iQ!  ~-W,y  Ja(t~x)  _  2  (t-x) 


/!» 


Q 


■  t  —  xj*  -r  \]  ’ 


from  (2.23)  we  obtain 


*/■-.!£ 


k(x,t ) 


y{  t)  dt  =  ~~  ,  -  a  <  x  <  a 


k(x,t)  =  1  J 


2\  +  n’(0)M0)  1  or !  A  ia[t_x) 


ha 


ia 


da  , 


(2.26) 


(2.27) 


(2.28) 


where  the  kernel  k(x,t)  is  bounded  in  the  closed  interval  —  a  <  (x,  t)  <  a  . 

From  (2.3b)  and  (2.21)  it  follows  that  (2.27)  must  be  solved  under  the  following  single¬ 
valuedness  condition: 


/  g{t)dt  =  0. 
J  -a 


(2.29) 


3.  Method  of  Embedding 

Even  though  the  direct  method  used  in  the  previous  section  yields  only  a  limited 
choice  for  representing  the  shear  modulus  y(y),  keeping  in  mind  that  the  mixed 
bounda-”  value  problems  which  are  of  practical  interest  tend  to  be  generally  local 
perturbeuion  problems,  it  may  not  be  difficult  to  approximate  f.i(y)  with  sufficient 
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accuracy  in  the  region  of  interest  by  using  the  representations  (2.10),  (2.11),  (2.12)  or 
(2.16).  To  solve  the  problem  one  could  also  follow  the  technique  developed  by  Varley 
and  Seymour  [10]  for  transforming  partial  differential  equations  with  variable 
coefficients  into  equations  with  constant  coefficients.  The  technique  is  more  general  in 
scope  but  also  is  more  complicated.  Following  [10]  we  first  replace  the  equilibrium 
equation  (2.2)  by 


M  (y)  Sx  "  dy  '  fL{-y' 


Introducing  the  functions  1  Y(x.y)  and  F(x.y)  which  satisfy 

d]V_d£  air  _  _dF  ,, 

dx  dy  '  dy  Ox  l  •  1 

we  now  express  w  and  /  in  terms  of  Il’(x,i/)  ,  F{x,y)  and  a  system  of  coefficient 
functions  wn  (y)  and  fn  ( y )  as  follows: 

Xu  f ).V  -  n 

w(xjj)  =  >  u'..tyi  - — v 11’, 

U  d’/~n 


Xu  O.V  -  n 

f(Fdj)  =  £/ n(i/)  Uf  V  -  n  f- 

u  ry 


Substituting  from  (3.3)  into  (3.1)  and  using  (3.2)  we  obtain 
a,/ 1 c hj)'  x-1" dyy-«{ dy ) ~ a/ 1 ay ) 


E  +  +  iv, 


(3.3a, b) 


av  /aw 

^  0  p,  v  l  ax 


,  uf  ds~n  (am  -f  £L(atf 

"atrv-nUx J~/oavvUx 


(3.4a,b) 
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An  admissible  set  of  functions.  //.  /„  and  wn  satisfying  ^3.4)  may  be  obtained  from 

fo  —  M  a'o  , 

In  ~  V  W*  ~  fl  “•*  -  I  '  71  =  1 . -V  • 

U/V  =  0. 

H  Wn  =  /n  -  f'n  -  1  '  77  =  1 .  -V  • 

/.V  =  o  .  ,3.5) 

From  3.3  it  may  be  shown  that 

,ro  '  'J:  -  •  f/’ji  =  v;'  Jt  ■ 

\  r1  '  ■// 

a’v  =  A  v  •  /v  =  F.v  •  (3-6) 


where  A"  v  and  I  v  are  constant. 

Equation  (3.2)  indicates  *  that  IT  is  a  harmonic  funtion  which,  for  y>0  . 
-  oc  <  x  <  oc  may  be  expressed  as 

J  ""  -4(a)  e~  'a  '  -J  ~ia"  da  .  (3.7) 

From  •  3.3)  and  '3.7)  it  follows  that 


±  [x 

2-J  -x 


: a  ;  1  ~ :a~  da 


where 

=  £  t»„(y)  (  -  !  a  i  )A  ",l  • 

0 

If  we  again  define  g(t)  by  (2.21),  from  (2.3)  and  (3.S)  we  obtain 


(3.3) 


(3.9) 


-lim  »■  !/.'  /"J  x  /Oh 

y-o.  —  y.a 


1/Y  ■  ■  c-  ■  a  !  y  +  .-Q(f-x) 


Q  /T 


;  a  Hi  0 .  a  | 


-  (7y.(.r.O)  =  a0(.n  .  -  a  <  x  <  a  . 


(3.10) 


By  separating  the  singular  kernel,  from  >3.10;  and  (3.9)  it  may  be  shown  that 


lim 

0 


i/:„ 


hix.'j.t)  at  =  — — 
-  hiQi 


a  <  x  <  a  , 


(3.11) 


hiz.y.t i  = 


■Wj  _ t_ 

u’.jiU)  if-.- 


m-  /"'*-**> 


,  -  !  a  j  i/  +  j'q(<  -  x) 


•'./  V!  a 


;a  //■  0 .  a 


(3.12) 


Bv  taking  the  limit  (3.12 1  becomes 


1  /”  J  f  1  1  (T  ■'  Z 1 

■=  /  r-1-;  -  £ci  j.  t '  f  o  -)  *  =  •  ~  a<  x<  a  , 

J  -  a  ^  t  £  ■'  ui  0  1 


(3.13) 


where  the  kernel  kc  is  bounded  in  -  a  <  [z.  t)  <  a  and  is  given  by 


k.'Z.  f'  = 


y  <  t°)  i  - 


1  „  ,.V  _  n 


Ja(t  -  x) 


V  wn[0)  (  -  jo'.  )A  ~  n 


£  <;0)  (  -  ■ 


V  at.  1  0)  (  -  n 


iin  q( f  —  x)  da 


(3.14) 


Equation  (3.13)  is  again  subject  to  the  condition  (2.29). 

As  an  example  consider  the  case  of  .V  =  1.  From  (3.5)  and  (3.6)  it  may  be  seen 


9- 


a’i  —  Kx  .  /[  —  £j  . 


(3.15) 


a/o  _  t'  ri  r 

(it/  -  A1  ^0  -  ^ 


(3.16) 


where  K{  and  £;  are  arbitrary  constants.  The  solution  of  (3.16)  may  be  obtained  as 


A',£,  >  0  : 


\T  ;'" 


A-  -\-r 


A',  £,  <  0  : 


(3.17) 


^  \  .  A'i  '  £,'  ;  1  y  ~  :/l ' 


L  ~  \j!  k  ‘‘"{£a  1’ 


(3.18) 


A'i  --  0  :  /j  =  -  I,  i  j  t  y;  j 


(3.19) 


1,^0:  /0=  - 

A  i 1  ;/  •*-  l/i 


(3.20) 


From  (3.6)  observing  that  /0  --  ^  ,  (3.17)  -  (3.20)  now  provides  four  different  functional 
terms  tor  wtiich  the  mixeu  boundary  value  problem  under  consideration  would  be 
analytically  tractable.  For  example,  assuming  Lx=  -  3^/tanhS,  Kx  =  -  3  tanhS 
and  Jf/j  =  £,  front  (3.17b)  we  find 


fij{y)  =  ^(y)  =  tanhiJy  —  d)/tann  6 


(3.21) 


VUI  y)  Ntzc 


iji/t  o  i/  cor 


(3.22) 


Also,  the  Fredholm  kernel  given  by  (3.14!  becomes 


kc(x.t)  =  2 3  ( coth2(6 )  -  1]  J 


sin  ait  —  x) 

3  +  a  coth{b)  Q 


(3.23) 


It  is  clear  that  solutions  for  /i(y)  with  greater  number  of  arbitrary  constants  may 
be  obtained  from  (3.5)  by  considering  cases  for  which  X  >  1  [10].  Technically  the 
greater  number  of  arbitrary  constants  would  provide  a  better  representation  of  the  given 
shear  modulus. 


4.  Examples 

We  first  consider  the  application  of  the  method  described  in  Section  2  and  give  two 
examples  where  we  assume  that  c0  >  0  or  cQ  —  0  (see  equations  (2.9)  -  (2.16)).  Thus, 
we  let  the  shear  modulus  have  the  form 

'  a)  n  [y)  =  //0  t1'1  ,  y  >  0.  (4.1) 


( b)  y(y)  =  y^(3y-  6j2  ,  y  >  0 


(4.2) 


where  3  may  be  positive  or  negative.  For  the  two  cases  (4.1)  and  (4.2)  the  Fredholm 
kernel  expressed  by  (2.2S)  becomes 


kj  i,fi 


-3 

o 


!  t  —  x 
t  —  £ 


1  1 

—  .r 


sm  ait—x) 

,J  a  (a  -  Ja*  -r  31  j 


da  , 


(4.3) 


t.f  i  =  - 

e 


d  f  x  sin  ait  -  x ) 


a 


da 


-  3  j  t  -  x  | 
2  6  t  -  x 


(4.4) 


From  (2.27)  and  (2. 28)  it  may  be  observed  that  the  dominant  part  of  the  integral 
equation  has  only  a  Cauchy  singularity.  Consequently  the  solution  of  (2.27)  has  the 
form  [11] 

g{t)  =  ,  -  a  <  t  <  a  ,  (4.5) 

n  a  -  tr 

where  G'  is  unknown  which  is  bounded  m  -  a  <  t  <  a  and  non-zero  at  t  —  ^  a. 


The  integral  equation  (2.27)  is  solved  numerically  by  normalizing  the  interval 


and  by  observing  that  its  fundamental  function  is  the  weight  of  Chebyshev  polynomials 
Tn(s).  Thus,  defining 


.  h(s)  .  G{as) 

t  =  as,  x  =  ar,  g(t)  =  o[s)  =  ,  ■  ,  h(s)  =  — —  , 

'l  1  —  s' 

The  unknown  function  !>{s)  may  be  expressed  as 


^  =  ^0)  ^ 


(4.6) 


(4.7) 


where  .40,  Aj,  .  .  .  are  unknown  constants  and  only  the  symmetric  loading  is  considered, 
that  is,  cr0  (x)  =  <70(  -  x)  and  w{x,y)  =  w( -x.y).  By  using  the  following  properties  of 
Chebyshev  polynomials 


0,  n  =  0,  -  1  <  r  <  1, 


I  f[  T"  (*)  ds 

J  -  1  ( s  —  1  -  s2 


<  C.'11.,(r).  n  =  1,2,...,  —  l<r<l, 

_  -==L^=  [•'-(!  r\/r)  'ir-l  ]  ,  n  =  0,1, 
•i  r  -  1  L 


(4.8) 


M  >i 


and  by  substituting  from  (4.6)  and  (4.7)  into  (2.27)  we  obtain 


E  -Ij  Pin 

-i(r)j  =  <?o(ar) 

,  -  1  <  r  <  1  , 

(4.9) 

0 

-  iM  =  !/  (  j  “  | 

ku  {ar,  as)  1 
kb  {ar,  as)  J 

ds  . 

'Jl-s2 

(4.10) 

Also,  by  using  the  orthogonality  conditions  of  Tn  (s)  from  (2.29),  (4.6)  and  (4.7)  it 
follows  that 


A0  =  0  ,  (4.11) 

The  simplest  way  of  solving  the  functional  equation  (4.9)  without  sacrificing 
accuracy  appears  to  be  truncating  the  series  at  n  =  N  and  using  an  appropriate 
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collocation  technique.  In  this  case  because  of  the  nature  of  the  problem  it  is  necessary 
to  increase  the  density  of  the  collocation  points  near  the  singular  points  r  =  1.  Thus, 

these  points  may  be  selected  as  follows: 


r.v  (rj  =  0  ,  r,  =  cos[  ^  (2i  -  1)  ,  i  =  1,  N  . 


(4.12) 


The  linear  algebraic  equations  obtained  from  (4.9)  by  substituting  r=ri;  i  =  l,  N 

would  then  give  the  unknown  coefficients  .dj,  .  .  .  ,  .dy. 

One  may  note  that  in  case  { b )  the  kernel  kb  is  known  in  closed  form  and  by  using 
(4.4)  and  (4.10)  the  functional  equation  (4.9)  may  further  be  simplified  as  follows: 


Y.  An  U2n  _2  (r)  +  JtYi  3*n(-n_  l)^i  =  ao(r)'  -  l<r<l 

where 

\  =z  ~  ,  r  =  cos9  . 
d 


(4.13) 


(4.14) 


Even  though  the  collocation  technique  mentioned  above  can  also  be  used  to  solve  (4.13), 
in  this  case  perhaps  a  more  appropriate  technique  would  be  to  use  a  weighted  residual 
method  and  reduce  (4.13)  to  an  infinite  algebraic  system.  Thus,  if  we  select  the  weight 
functions  U2n  (r)  'll  —  r2  and  use  the  orthogonality  conditions 


l  1  i - j  ^ '  !  ^  j 

A  U,  (r)  Uj  ( r)  'll  -  i-2  dr  = 

J  -  1  i  i  =  j 

\  **' 


(4.13)  may  be  reduced  to 


(4.15) 


■4,  -  £  C<J  -4;  +  di-i  ,i=  1,2,  .  . 


(4.16) 


where 


d,  =  -  J  J[{  <?q  (r)  Uu  (r)  'J  1  -  r1  dr  ,  i  =  0,1,2, 


(4-17) 


(2i  —  2)2  —  (2j  —  1  )2  (2i)2  —  (2j  —  1) 


(4.18) 
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The  infinite  system  (4.16)  can  be  solved  by  reduction  and  convergence  is  assured 
provided  the  system  is  regular  [12] .  The  condition  for  regularity  is 


£  I  ctJ  |  <1,  i  —  1.2,  .  .  . 

J  =  1 

Performing  the  infinite  sum,  (4.19)  may  be  expressed  as 

2Ai*  ~ H  i  9  y''  ( _ 1 _ 1 _ 

(4t-3)(4i-l)  ‘;t'1U2»-2^-(2i-l)2  (2i)2  —  (‘2j  —  l)2 


(4.19) 


<  1. 


(4.20) 


In  (4.20)  the  bracketed  quantity  is  greatest  for  1=1.  Hence  the  condition  of  regularity 
becomes 

A  <  if  or  T  <  if  '  (4-21) 

For  A  >  3-/16  the  system  is  quasi-regular  in  the  sense  that  there  exist  an  i  for  which  the 
infinite  system  beginning  with  the  (i  + 1)  st  equation  is  regular.1 

In  the  problem  under  consideration  the  physically  important  quantities  are  the 
crack  opening  displacement  udx.O)  and  the  shear  stress  <ry.(x, 0)  along  the  interface, 
particularly  the  stress  intensity  factor.  From  (2.21),  (4.6)  and  (4.7)  ut(x,0)  may  easily 
be  obtained  as  follows: 

u'(x.O)  =  [  g{t)  dt=  --hr,  Y.  u  -I..  L'  in-*  (  1)  ^  -  **  ,  M  <  a  .  (4.22) 

J  -a  Mu  J  y 

Referring  to  (2.23)  and  (2.27)  it  may  be  observed  that  the  left  hand  side  of  the  integral 
equation  represents  <7y.(x,0)  for  |x|  >  a  as  well  as  for  |x|  <  a.  Thus,  by  substituting 
from  (4.6)  -  (4.S)  into  (2.27),  for  x  =  ar  and  |r  \  >  1  we  obtain 

^here  seems  to  be  a  curious  analogy  between  the  antiplane  shear  problem  for  a  half  plane 
with  the  shear  modulus  given  by  (4.2)  and  the  cover  plate  problem  described  in  [13],  In  [13]  the  cover 
plate  was  approximated  by  a  membrane  of  thickness  h  and  elastic  constants  /q  and  tq  ,  and  A  was  the 
stiffness  parameter  defined  by  A  =  /q  ( 1  +  fq)/[-Mih  (l  +  «2)]  where  n2  and  k2  are  the  elastic 
constants  of  the  substrate.  In  the  two  problems  the  analogous  physical  quantities  are  the  interface 
shear  stress  crjy  (in  [13])  and  p.  (0)  0  (s)  (sec  (2.21)  and  (4.G). 
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Based  on  simple  physical  considerations  it  may  be  observed  that  the  solution  given 
in  this  study  for  a  nouhomogeneous  half  plane  having  a  shear  modulus  n(y)  ,  y> 0,  is 
valid  also  for  the  corresponding  infinite  medium  in  which  y  =  0  is  a  plane  of  material 
symmetry,  that  is.  /h  -  y)  =  /by),  -  rc  <  y  <  tx:.  Thus,  in  the  examples  considered  if  we 
let  the  material  nonhomogeneity  parameter  3  be  zero  in  (4.1),  (4.2)  and  (4.26),  n 
becomes  a  constant  /n,0)  and  the  problem  reduces  to  that  of  a  homogeneous  plane 
having  a  crack  which  is  subjected  to  symmetric  crack  surface  tractions.  In  fact  from 
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(2.17),  (2.27)  and  (2.2S)  it  may  be  seen  that  for  3  —  0  we  have  A  =  |  a  |  and  /i'(0)  =  0 
giving  k(x,t)  =0.  For  example,  if  <ru(x)  =  -  pu,  from  (4.6) -(4.11)  and  (4.23)  it  follows 
that, 


M2n  -  i  —  0,  •41 


Po <  An  =  0.  in  >  1),  g{x) 


<ryt(x,0) 


—  Po  'ia  ■ 


Po 

Po 


X 

\a‘  -  x2 


;  x  |  <  a, 


(5.1) 


Similarly,  observing  that  r  =  4-  f'0( r)j/4  ,  for  a  parabolic  traction 

<t0(  j)  =  -  p2{x/a)2  it  may  easily  be  seen  that  =  A2  =  p2  /4  and 


cry3(  j.U) 


Pi  A 
N  x2  -  a2 


X3] 

)  |  ,  x  >  a  , 


(5.2) 


giving  ki  =  p2\a/ 2  (see  Table  1).  Also,  in  the  third  example  considered  in  Section  3, 
from  (3.23)  it  is  seen  that  for  3  =  0  kcLv.t)  =  0  and  (3.13)  reduces  to  an  integral 
equation  for  a  homogeneous  material  having  the  shear  modulus  /z(0). 

For  a  plane  with  various  forms  of  nonhomogeneities  considered  in  Sections  2  and 
3  the  calculated  stress  intensity  factors  are  given  in  Tables  2  and  3.  In  these  results  pQ 
and  p2  again  refer  to  the  loading  parameters  defined  in  Table  1.  In  a  limited  way  the 
results  for  uniform  loading  arc  al-o  displayed  in  Figures  2-4.  The  materials  considered 
in  these  three  examples  exhibit  some  distinct  features.  The  hyperbolic  functions  used  in 
Table  2  and  Fig.  2  to  represent  ..e  shear  modulus  have  constant  asymptotic  values  for 
y-*7c,  namely  p0,'tanh26  and  p0/ cot.h2S.  where  p0  =  /z(0)  (see  the  insert  in  Fig.  2).  In  the 
parabolic  distribution  for  /<(?/)  assumed  in  Fig.  3  and  Table  3,  for  3  <  Op  becomes  zero 
at  y  =  6 / 3  and  hence  the  results  correspond  to  a  nonhomogeneous  layer  of  finite 
thickness  with  vanishingly  decreasing  shear  modulus  away  from  the  interface.  For  this 
reason  in  this  example  the  stress  intensity  factors  corresponding  to  /?  <  0  are 
considerably  greater  than  that  of  the  other  examples  shown  in  Figures  2  and  4  and 
Tables  2  and  3  in  which  u[y)  is  also  a  monotonicallv  decreasing  function. 

Let  us  now  examine  the  asymptotic  behavior  of  the  stress  intensity  factor  as 
3-*  —  x.  First,  from  (4.6).  (4.7)  and  (4.25)  it  can  be  seen  that  k3(a )  may  also  be 
expressed  in  the  following  alternate  form: 
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ki(a)=-Uma  //(O)  |2(u  -  x)  g(x)  . 


(5.3) 


In  the  examples  given  in  Table  3  and  Figures  3  and  4,  it  may  now  be  seen  that  as 
+  o c  the  stiffness  of  the  medium  would  increase  indefinitely  and.  under  loadings  of 
finite  magnitude,  the  crack  opening  displacement  or  g(x).  and  consequently,  k3(a)  would 
tend  to  zero,  these  physically  expected  trends  may  be  observed  in  Figures  3  and  4  and 
in  Table  3.  Similarly  in  the  example  given  in  Table  3  and  Fig. 4  for  y(y)  =  fiQexp(23y), 
the  stiffness  of  the  medium  would  decrease  indefinitely  as  3-*  -oo,  and  consequently 
k3(a)  would  tend  to  infinbv.  Also,  for  3  <  0  and  p(y)/p0  =  (3y  +  S)2/62  the  “thickness” 
of  the  medium  yQ  =  —  Si 3  and  as  a  result  its  stiffness  would  decrease  indefinitely  as 
3-*-  x,  again  causing  k3'  a)  to  approach  infinity. 

The  examples  considered  in  Table  2  and  Fig.  2  are  quite  different  and  somewhat 
more  practical  than  the  examples  shown  in  Table  3  in  that  for  /3-» oo  the  medium 
becomes  homegeneous  with  a  shear  modulus  /<(0) jcoth2S  or  p(0)/tanh2S.  Thus,  in  the 
limiting  case  of  3  —  x  the  stress  intensity  factor  ought  to  be  independent  of  the  shear 
modulus  and  should  have  the  values  given  in  Table  1.  This  trend,  of  course,  is  not 
observed  in  Table  2  arid  Fig.  2.  The  explanation  of  the  discrepancy  may  be  found  in 
the  fact  that,  becouse  of  the  discontinuous  change  taking  place  in  /x(0)  at  f3  =  o o,  the 
only  physical  quantity  that  should  be  expected  to  vary  continuously  as  j3-+oo  is  the 
strain  energy  release  rate  which  may  lie  defined  by  and  expressed  as 

-  «'«•»)<»*.  (*•*> 

where  E0  is  the  externally  added  or  internal}'  released  energy  that  is  available  for 
fracture.  Physically  the  problem  may  be  visualized  as  a  crack  problem  for  a 
homogeneous  half  plane  bounded  to  a  rigid  substrate  through  a  thin  nonhomogeneous 
layer.  The  objective  is  to  determine  the  asymptotic  behavior  of  the  solution  as  the  layer 
thickness  approaches  zero. 

To  gain  a  better  insight  to  the  problem  we  consider  the  simpler  example  of 
piecewise  homogeneous  medium  shown  in  Fig.  5.  The  mode  III  crack  problem  described 
by  Fig.  5  is  that  of  two  idendical  homogeneous  half  spaces  having  the  shear  modulus 
which  are  bonded  through  a  homogeneous  layer  of  shear  modulus  and  thickness  2 h. 
The  layer  contains  a  crack  of  length  2 a  in  the  plane  of  the  symmetry.  Following  a 
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technique  described  in  section  2.  the  mixed  boundary  value  problem  may  easily  be 
reduced  to  the  following  singular  integral  equation  [14]: 


i  I  _a{  l~x  ~  p  *  k^x<  *)}  ^  dt  ~  ^FT  ’  ~a<x<  a  < 

fj.  *  =  -  -j-K  k\(x,t)  =  2  f  ' - - — —7 - sina(t  —  x)da 

Pi  t/V  ’  J  0  /i  ‘  e  ” 2a '  +  1 


— 2a/i - sina(t  —  x)da 


(5.5a,b) 


First,  we  simply  note  that  for  h  =  oo  the  problem  becomes  one  of  a  cracked 
homogeneous  plane  with  shear  modulus  /q  for  which  the  solution  is  given  by  (5.1)  and 
(5.2)  and  Table  1,  specifically 


_  ,  Ma)  ,  Q(«)  ,  „  trpoa 

<7°(-r)  ~  P°:  ~L  §Q  ~  h  5o  “  4/q  ' 

a  lx)~  -  D(I)2.  Mul  - 1  M_0o3  q  -IE?" 

J  ~  p2',7i  “2'  “  U-  ’  yi  “  4/q  ’ 


(5.6a,b) 


Referring  to  (5.4)  we  note  that  the  strain  energy  release  rate  Q  calculated  in  this  article 
and  given  in  various  figures  and  tables  and  in  equations  (5.4),  (5.6)  and  elswhere 
correspond  to  one  quarter  of  the  total  enwrgy  released  in  a  cracked  medium  such  as  that 
shown  in  Fig. 5.  In  the  other  limiting  case  of  /t~» 0  it  may  easily  be  shown  that 


and  (5.5a)  becomes 

I  r  Mdt  =  .  -  a 

"  J  -a  f-  2  P2 


<  x  <  a  , 


which  is  the  integral  equation  for  a  homogeneous  medium  having  the  shear  modulus  /q. 
Thus,  from  (5.S)  it  follows  that 


,rU  n  .  £f0  _  ,  r>  _  _  P±ft 

<To(X)  -  ~Po-  §2  ~1'  ~  4/q  “ 


/  x  /.rV  YW  no-  p  TP2a  Pip 

a°(x)  =  -  P2{a )  ■  =  °-2^  §3  =  ^  =  /qQi- 


(5.9a, b) 
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Referring  now  to  Fig. 5,  we  observe  that  as  long  as  h  >  0,  the  stress  intensity 
factor  and  the  strain  energy  release  rate  are  defined  by 


k3(a)=-lirna  p(0)  ]2(a  -  x)  Q(a) 


(5.10) 


where  g(x)  is  obtained  from  (5.5a).  For  the  loading  conditions  given  in  Table  1  and  for 
Pi/P2  =  2  and  fil/[i2  =  0.o  these  calculated  results  are  given  in  Table  4  and  are  also 
partially  shown  in  Fig. 6.  The  limiting  values  of  Q(a)  for  h-*o c  and  h-*Q  and  fc3(a)  for 
h-*z c  are  clear  and  unambiguous  and  are  given  by  (5.6)  and  (5.9),  respectively.  From 
(5.10)  (5.9)  and  (5.6)  it  is  also  seen  that  for  h  =  +0 


~k\{a) 

4/ti 


9(a)  = 

Q(a)  = 


rp-u 


1  '  t-l  1  f- 

4p,  •  =\jr,Pu'r< 

=  0.25  =  =  0.25  .  t3(a )  =  J  ^jp2Na. 


(5.11a,b) 


Table  4  and  Fig. 6  show  that  these  arc  indeed  the  limiting  values  of  the  results  obtained 
from  (5.5a)  and  (5.10). 

Considering  the  fact  that  for  /;-* 0  the  problem  reduces  to  that  of  a  homogeneous 
plane,  the  limiting  value  of  b3  still  appears  to  be  somewhat  paradoxical.  The 
explanation  of  the  apparent  discontinuity  in  k2  at  h—+  0  lies  in  the  fact  that  for  very 
small  k,  k3  obtained  from  the  integral  equation  (5.5a)  and  shown  in  Table  4  is  the 
measure  of  near  field  quantities.  Far  field  quantities  should  ,of  course,  be  that  of  the 
homogeneous  infinite  medium.  Without  a  lengthy  singular  perturbation  analysis  this 
can  be  shown  by  examining,  for  example,  the  crack  opening  displacement  u;(x,0)  very 
near  the  crack  tip  as  h  approaches  zero.  The  results  calculated  from  (5.5a),  (4.5)  and 
(4.22)  are  shown  in  Figures  7  and  S.  Figures  7a  and  Sa  show  u;(x,0)  and  (?(x)  for  various 
values  of  p l/;z2  and  hja.  From  Fig.  7a  it  is  seen  that  away  from  the  crack  tip  as  /i-»0 
tr(x,0)  converges  to  the  corresponding  homogeneous  infinite  plane  values  (indicated  by 
h  =  0).  Fig.  7b  shows  a  magnified  version  of  u;(x,0)  near  the  crack  tip  for  —  2. 

Note  that  tc  ~  p0/p.  The  curves  marked  by  ( hja  =  0,  p2/Mi  =2)  and  by  fil  =  ji2  are 
simple  ellipses.  As  expected,  away  from  the  crack  tip  the  curves  marked  by  small 
nonzero  values  of  hja  are  seen  to  converge  to  ///a  =  0  ellipse.  However,  near  the  crack 
tip  the  behavior  of  these  curves  is  quite  different  and  as  h-* 0  they  seem  to  have  a 
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different  slope  than  that  of  the  corresponding  ellipse.  These  may  be  seen  somewhat 
better  in  Fig.  3  showing  G{x),  G(a)  being  the  measure  of  k3(a)  (see  Eq.  (4.5)).  From 
(5.10a)  note  that  for  h  =  0,  that  is.  for  the  homogeneous  plane  /r  =  ji2, 
k3  =  -n2G{a)/-{a,  whereas  for  h  >  0  k3  =  —  filG(a)/\a.  Hence,  the  straight  line 
corresponding  to  h  =  0  gives 


G{x) 


1  Pox 

2  Pi 


P  of 
fl2  - 


=  A)N^- 


(5.12) 


The  values  of  G(a)  for  nonzero  values  of  h  shown  in  the  figure  give  the  stress  intensity 
factors  listed  in  Table  4.  The  figure  snows  that  at  h  =  +0  G(x)  would  have  a  spike  the 
value  of  which  seems  to  approach  the  independently  calculated  value  of  \nl/ n2.  The 
distribution  of  G{x )  also  shows  the  expected  “pinching"  of  the  crack  opening 
displacement  very  near  the  crack  tip,  as  exhibited  by  two  inflexion  points. 

Figures  9  and  10  give  some  sample  results  showing  the  influence  of  the  material 
nonhomogeneity  parameter  3  on  the  crack  opening  displacement  u,'(x,0),  -a<x<a.  As 
3  increases,  in  both  examples  the  material  stiffness  also  increases.  Consequently,  in 
both  cases  the  crack  opening  displacement  is  seen  to  decrease  for  increasing  values  of  i3. 

Based  on  the  results  given  in  Tables  2  and  3  one  may  conclude  that  in  fracture 
mechanics  analysis  involving  nonhomogeiieous  materials  the  “crack  driving  force"  (as 
measured  by  the  stress  intensity  factors  or  the  strain  energy  release  rate)  may  differ 
quite  considerably  from  the  values  given  by  the  corresponding  homogeneous  materials. 
One  may  also  note  that  in  diffusion  problems  the  square-root  nature  of  the  flux 
singularity  at  the  crack  tip  is  not  influenced  by  the  material  nonhomogeneity. 
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Table  1.  Loading  conditions  used  and  the  corresponding 
stress  intensity  factors  for  a.3— 0. 


m 

~Po 

-  p2(x/a)~ 

k'o(a) 

| 

Port 

ljP2rt 

Table  2.  The  variation  of  stress  intensity  factors  for  various 
loading  conditions  shown  in  Table  1.,  8  =  1. 


p(y\'p  o 

coth~(3y  +  8 ) 
co th'(6) 

p(y)/p o  = 

C) 

tanh~(3y  +  6 ) 
tanh2{8) 

3a  i 

i 

i 

*3 

Po\a 

k-i 

Pi^a 

\ 
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Table  3.  The  variation  of  stress  intensity  factors  for  various 
loading  conditions  shown  in  Table  1.,  <5  =  1. 
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Table  4.  The  stress  intensity  factors  and  the  strain  energy 
release  rates  in  two  homogeneous  '  f  planes  bonded 
through  a  homogeneous  layer  containing  a  crack  (Fig.  5). 
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Fig. 2.  Stress  intensity  factor  and  Strain  energy  release  rate  for 
a  nonhomogeneous  half  plane  under  uniform  antiplane  shear 
loading,  cryi(x,0)  =  -  p0,  5o  =  ^Po0/4/-^, 

A:  n(y)  =  n0tanh2(Py  + 8)/tanh26, 

B:  p(y)  =  n0coth2(/3y  +  6 )/coth26,  6  =  1. 
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Fig. 3  Stress  intensity  factor  for  a  nonhomogeneous  half  plane 
subjected  to  uniform  crack  surface  traction 

o't/z^-O)  =  -p0,  shear  modulus  p(y)/p0  =  (£y  +  <5)2/<52,  <5  =  1. 


Fig. 5  The  symmetric  modem  crack  problem  for  two 
homogeneous  half  spaces  bonded  through  a  homogeneous  elastic 
layer. 


h/a 


Fig.  6  Stress  intensity  factor  and  strain  energy  release  rate 
as  a  function  of  layer  thickness  to  crack  length  ratio  , 
§o  =  *p%a/ Apx, 
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Fig. 7a  Crack  surface  displacement  in  a  layer  bonded  to  a 


semi- infinite  homogeneous  plane  under  uniform  crack 
surface  shear  loading  ayZ(x,Q)  —  —  p0  ,  wQ  =  - . 


Fig. 7b.  Crack  surface  displacement  in  a  layer  bonded  to  a 
semi-infinite  homogeneous  plane  under  uniform  crack 
surface  shear  loading  <TyZ(x,0)  =  —  p0  ,wQ  =  ^-p—. 
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Fi0  )  Crack  surface  displacement  in  a  semi-infinite 
nonhomogeneous  plane  under  uniform  crack  surface  shear 
loading  cryi(x,0)  =  -  p0;  fi{y)/fi0  =  tanh2{py  +  6)/tanh2(6 ), 
S  =  l,u>0  =  apJfiQ. 
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Fig.  10  Crack  surface  displacement  in  a  semi-infinite 
nonhomogeneous  plane  under  uniform  crack  surface  shear 
loading  <xy2(x,0)  =  -  p0;  /i(y)/p0  =  ,w0  =  ap0/p0. 
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